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Abstract 

Vector boson plus jet production is interesting for Higgs search, beyond the Standard Model physics and provides 
standard candles for calibration. This is complementary to inclusive jet production measurements which provide 
precision tests of perturbative QCD. A multitude of W/Z plus heavy and light flavour jet measurements in pp collisions 
at a centre of mass energy of -\fs = 1 .96 Te V is discussed. Next-to-Leading order perturbative QCD predictions and 
various models are compared to the measurements. 



1. Introduction 

The discussed measurements are accomplished by the 
multi-purpose collider experiments D0 and CDF with 
their broad particle identification capabilities through 
central tracking, fine granulated calorimeters and muon 
spectrometers. The detectors are located at the Tevatron 
accelerator at Fermilab, where proton anti-proton colli- 
sions take place at a centre of mass energy of 1 .96 TeV. 
Vector boson plus jet production at the Tevatron is com- 
plementary to the kinematic regimes of the HERA ac- 
celerator and fixed target experiments. It provides han- 
dles to Parton Distribution Functions (PDF), initial and 
final state gluon radiation and the transverse momen- 
tum (pr) spectra of the jets and the vector bosons. The 
measurements provide standard candles for calibration 
and tuning of Monte-Carlo (MC) event generators. The 
understanding of these processes is important for Stan- 
dard Model (SM) as well as Beyond the SM (BSM) phe- 
nomenology. The Tevatron dataset is now large enough 
to confront predictions and it has a unique kinematic 
overlap with the LHC and the expected SM Higgs mass 
range. 

The data of the measurements are fully corrected for 
instrumental effects. Therefore they can be directly used 
for testing and improving existing MC event genera- 
tors and any future calculation and model. Next-to- 
Leading Order (NLO) and Leading Order (LO) per- 
turbative QCD (pQCD) predictions from MCFM lH 
are compared to data taking non-perturbative effects 
of hadronisation and underlying event from simulation 



into account in the prediction, i.e. the comparison takes 
place at the hadronic final state or synonymous at the 
particle level. The relative uncertainties of the mea- 
surements are dominated by the Jet Energy Scale (JES) 
uncertainty. Iterative seed-based infrared safe midpoint 
cone jet algorithms are used by the D0 [2] and CDF 
experiments in Run II. 

In the following sections 2-6 the Tevatron vector bo- 
son plus light flavour jet measurements are presented, 
followed by the vector boson plus heavy flavour jet mea- 
surements in sections 7-11. 



2. D0 inclusive Z/y*(-» ee) + jet cross section 

The Z + jet production cross section is measured J3l 
differentially in inclusive pj bins of leading, second and 
third leading jet, making use of an integrated luminosity 
of l.Oftr 1 . Electron/positron candidates above a trans- 
verse momentum of pj > 25 GeV are selected in an 
absolute pseudorapidity range of |77 e | < 1.1, 1.5 < \rj e \ < 
2.5 to form Z boson candidates in a mass window of 
65 < M (e < 1 15 GeV. Jets with a cone radius of R = 0.5 
and a prf > 20 GeV are selected in a rapidity range of 
|vj e t| < 2.5. The cross section is normalised to inclusive 
Z + jet production which cancels uncertainties on lu- 
minosity and most of electron trigger and identification 
uncertainties. The phase space of the selected events is 
extrapolated to the full lepton kinematics. While the 
(N)LO pQCD predictions are in agreement with data 
there are large differences between the different mod- 
els of PYTHIA |5|. HERWIG [6] (+ JIMMY), ALP- 
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GEN 0] + PYTHIA and SHERPA ft The experimen- 
tal errors are small and dominated by statistics, allowing 
for future improvements. 



3. CDF inclusive Z/y*(-» ee) + jet cross section 

The Z + jet production cross section is measured |@] 
differentially in inclusive transverse momentum bins of 
the leading and second leading jet, making use of an in- 
tegrated luminosity of 1.7 fbr 1 . Electron/positron can- 
didates above a transverse energy of E e T > 25 GeV are 
selected in an absolute pseudorapidity range of \rj e \ < 
1.0, 1.2 < \r} e \ < 2.8 to form Z boson candidates in a 
mass window of 66 < M ee < 1 16 GeV. Jets with a cone 
radius of R — 0.7 and a p^ > 30 GeV are selected in a 
rapidity range of |yj et | < 2.1. There is good agreement 
between data and NLO pQCD prediction. 



4. D0 inclusive Z/y*(-» ////) + jet cross section 

Differential angular distributions of the Z + jet produc- 
tion cross section is measured 1 10] in inclusive bins of 
transverse Z boson momentum making use of an inte- 
grated luminosity of 1 .0 fbr 1 . Oppositely charged muon 
candidates above a threshold of pj > 15 GeV are se- 
lected in an absolute rapidity range of |yj < 1.7 to 
form Z boson candidates in a mass window of 65 < 
M m < 115 GeV. Jets with a cone radius of R = 0.5 
and a pf > 20 GeV are selected in a rapidity range 
°f lyjetl < 2.8. The differential cross section in two in- 
clusive bins of Z boson transverse momentum (pj > 
25 GeV and pj > 45 GeV) is measured as a function of 
the angular variables rapidity sum, rapidity difference 
and the azimuthal angle difference between the leading 
jet and the Z boson. ALPGEN 01 and SHERPA M 
include up to three partons in the matrix element cal- 
culations. The binning is chosen such, that the detector 
resolution causes little migrations between bins. There 
is less agreement between data and predictions in the 
lower pj > 25 bin. This holds for the non-pQCD as 
well as for the pQCD predictions. Among the non- 
pQCD predictions SHERPA (1.1.3) provides the best 
description of the angular distributions, in particular 
Ay(Z, j), in the p z T > 45 GeV bin. In the pj > 25 GeV 
bin the measured cross section is <x(Z + jet)/cr(Z) = 
[47 ± l(stat.) ± 2(syst.)] ■ 10~ 3 compared to the pQCD 
NLO prediction of [40 + 3(scale) ± l(PDF)] • 10~ 3 and 
the pQCD LO prediction of [40 ± 8(scale) ± l(PDF)] • 

io- 3 . 



5. CDF Z(-> ee,fifi) + 1 jet p T balance 

CDF measured the Z boson + 1 jet exclusive pj bal- 
ance jlin . making use of an integrated luminosity of 
4.6 fbr 1 and sets precision limits on the measurements 
of Standard Model (SM) jet and vector boson ob- 
servables, relevant for the discovery potential of new 
physics. Oppositely charged lepton candidates {( = 
e,p) above a threshold of Pj,E T > 18 GeV are se- 
lected to form Z boson candidates in a mass window of 
80 < M m < 100 GeV. A leading jet with a pf > 8 GeV 

is selected in a pseudorapidity range of 0.2 < |77 det | < 0.8 

jet 

to avoid cracks in the central calorimeter. The measure- 
ment is accomplished for three different cone algorithm 
radii R = 0.4, 0.7, 1 .0. The azimuthal angle between the 
Z boson and the jet has to satisfy Acp > 3.0 rad to pro- 
vide back-to-back configurations in the plane transverse 
to the beam axis. Predictions for Z boson + jet produc- 
tion are obtained from PYTHIA and ALPGEN 10]. 
The CDF JES is calibrated via tuning of the calorimeter 
response to single particles. Therefore the Z boson plus 
jet pj balance provides an independent test of the CDF 
JES. While the measured track in jet distributions agree 
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Figure 1 : Leading jet plus Z boson pj balance data over simula- 
tion (PYTHIA) ratio as a function of the azimuthal angle difference 
between the leading and the second leading jet. Data has more sub- 
leading jets close in angle to the leading jet. 

with the quark and gluon jet fractions of PYTHIA 01, 
a higher rate of sub-leading jets collinear to the leading 
jet is observed in data (see Fig. [T]). This does set lim- 
itations on the precision of the JES determined at the 
ATLAS and CMS experiments by means of the jet plus 
Z boson pj balance of the order of 3%. Only large an- 
gle final state radiation observed in form of sub-leading 
jets turns out to be able to explain the discrepancy. 
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6. CDF inclusive W(-* ev) + n=l-4 jets cross section 



8. CDF inclusive cr(Z + 6)/tr(Z) ratio (Z -> 



The W + n jet inclusive production cross section (« = 
1 - 4) is measured lITill differentially in inclusive trans- 
verse momentum bins of the nth leading jet, making 
use of an integrated luminosity of 320 ptr 1 . Elec- 
tron/positron candidates above a transverse energy of 
E e T > 20 GeV are selected in an absolute pseudorapid- 
ity range of \rj e \ < 1.1. A missing transverse energy of 
E j > 30 GeV and a transverse mass of nij > 20 GeV is 
required to select W boson candidates. Jets with a cone 
radius of R = 0.4 and a p^ > 20 GeV are selected in a 
pseudorapidity range of |77j e tl < 2.0. The measurement 
benefits of a production cross section which is about ten 
times higher than the one for Z + jets production. At 
the same time the multi-jet and top quark production 
background needs to be controlled. There is good agree- 
ment between data and the NLO pQCD prediction. At 
low transverse jet momentum the MC event generators 
ALPGEN 0] plus PYTHIA H (with MLM matching 
of the matrix element to the parton shower) and MAD- 
GRAPH H! plus HERWIG J3] (with CKKW match- 
ing) need a better modelling of the underlying event. 



The inclusive Z + b jet production cross section is mea- 
sured [15] as the ratio over the inclusive Z production 
cross section, making use of an integrated luminosity of 
2.0 fbr 1 . Z boson decays into a pair of charge conju- 
gated electrons or muons are considered. The cross sec- 
tion is measured differentially in bins of jet and heavy 
flavour jet multiplicity, as functions of leading b jet 
transverse energy and pseudorapidity and as a func- 
tion of the Z boson transverse momentum, b jets are 
identified by means of a secondary vertex tagging al- 
gorithm. The leading and second leading leptons (e,p) 
are required to satisfy \/j ,Ej > 18 GeV and p'p, Ej > 
10 GeV. The invariant dilepton mass has to lie in the in- 
terval 76 < Met < 106 GeV. Jets with a cone radius of 
fl = 0.7anda£f>20GeV are selected in a pseudora- 
pidity range of |?7j et | < 1.5. The measured cross section 
ratio cr(Z+b)/o-(Z) = [3.32 ± 0.53(stat.) ± 0.42(syst.)] x 
10" 3 is in agreement with NLO pQCD predictions. Both 
data and theory have large uncertainties. While data un- 
certainties are statistics dominated the large uncertainty 
of the theory prediction arises from the missing NLO 
pQCD prediction for Z + bb production which in turn 
causes a large scale dependence. 



7. D0 inclusive cr(Z + b)/(r(Z + j) ratio (Z -> ee, fifi) 



The inclusive Z + b jet production cross section is mea- 
sured as the ratio over the inclusive Z + jet pro- 
duction cross section, making use of an integrated lu- 
minosity of 4.2 fbr 1 . Z boson decays into a pair of 
charge conjugated electrons or muons are considered. 
Muon candidates with > 10 GeV in the pseudora- 
pidity range of l^tl < 2.5 and electron/positron can- 
didates with p e T > 15 GeV in the pseudorapidity range 
of 1 77ci et I < 2.0 are selected together with jets (R = 0.5). 
The leading jet has to have a transverse momentum of 
p>* > 20 GeV in a pseudorapidity range of 1 77det I < 1-1. 
The b flavour jets are separated from c and light flavour 
jets by means of a neural network jet tagging algorithm 
based on the longer mean lifetime of heavy flavour 
hadrons with respect to light ones. The flavour frac- 
tions are determined by a fit of lifetime variable tem- 
plates to the data. The templates are obtained from 
ALPGEN H plus PYTHIA H while the cross sec- 
tions are taken from NLO pQCD calculations. The mea- 
sured total cross section ratio is cr(Z + b)/cr(Z + j) = 
0.0176 ± 0.0024(stat.) ± 0.0023(syst.) consistent with 
the NLO pQCD prediction of 0.0184 + 0.0022. 



9. D0 inclusive <r{W + c)/(r(W + j) ratio (W -> tv) 

The inclusive W + c jet production cross section is mea- 
sured lfl6ll as the ratio over the inclusive W + jet pro- 
duction cross section, making use of an integrated lumi- 
nosity of 1.0 fbr 1 . The leptonic W boson decay with 
a muon/electron in the final state is considered. The 
cross section is measured differentially in bins of jet pj. 
Lepton candidates with p c T > 20 GeV and events with a 
missing transverse energy of Ej > 20 GeV are selected. 
Jets with a cone radius of R — 0.5 and p 1 ^ 1 > 20 GeV are 
selected in a pseudorapidity range of |//j et | < 2.5. W bo- 
son + c quark production is sensitive to the s quark PDF 
content up to scales of Q 2 = 10 4 GeV 2 . The charge sign 
of the lepton from the W boson decay is opposite (OS) 
to the one of the c quark. This information is exploited 
in tagging an oppositely charged muon in the jet from 
the c quark fragmentation. The same sign (SS) sub- 
tracted background is small (~ 1%). The multi-jet back- 
ground is determined from data. The measured cross 
section is corrected for acceptances and efficiencies. It 
amounts to cr(pp — > W + c-}et)/cr(pp — > W + jets) = 
0.074 + 0.019(stat.)+° °} 2 (syst.) and is within the large 
errors in agreement with the prediction of 0.044 ± 0.003 
from ALPGEN 0] plus PYTHIA H. 
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10. CDF semi-incl. W + c jet cross section (W -* tv) 

The W + c jet exclusive production cross section is mea- 
sured [ 17] for selected events with one or two jets, mak- 
ing use of an integrated luminosity of 1 .8 fb 1 . Lep- 
ton candidates with p e T > 20 GeV in the pseudorapid- 
ity range of \rj(\ < 1.1 and events with a missing trans- 
verse energy of $Y > 25 GeV are selected. Jets with 
a cone radius of R — 0.4 and > 20 GeV are se- 
lected in a pseudorapidity range of \rjj et \ < 1.5. Charm 
flavour jets are tagged by means of a soft muon in the 
jet. The OS-SS subtracted data distributions are in good 
agreement with the prediction of ALPGEN 10] with 
PYTHIA fl. The total measured cross section amounts 
to oiW(->tv)+c) = 9.8±2.8(stat.)+|-4(syst.)+0.6(lum.)pb 
and is in good agreement with the NLO pQCD predic- 
tion of aiW(-my) +c) = 1 1.0+14 pb for pp et > 20 GeV 



and | j], 



c-jet I 



< 1.5. 



11. CDF semi-incl. W + b jet cross section (W -* tv) 

The W + b jet exclusive production cross section is mea- 
sured [18] for selected events with one or two jets, mak- 
ing use of an integrated luminosity of 1.9 fb -1 . The 
leptonic W boson decay with a muon/electron in the fi- 
nal state is considered. Lepton candidates with p e T > 
20 GeV and |%| < 1.1 and events with a missing trans- 
verse energy of tyj > 25 GeV are selected. Jets with 
a cone radius of R = 0.4 and pf > 20 GeV are se- 
lected in a pseudorapidity range of |?7j et | < 2.0. Bottom 
flavour jets are tagged by means of a secondary ver- 
tex tagging algorithm. A maximum likelihood fit to the 
vertex mass distribution of tagged jets is applied to ex- 
tract the flavour fractions of the selected data sample. 
The flavour templates are obtained from ALPGEN [0] 
with PYTHIA m and MADEVENT JH] for single top 
quark production. The measured total cross section is 
aiW^ty) +b) = 2.1 A ± 0.27(stat.) ± 0.42(syst.) pb. This 
has to be compared to the NLO pQCD prediction of 
cr(W { ^ {v) + b) = 1 .22 ± 0. 14 pb and the LO prediction of 
oiW^y) +b) = 0.78 pb from ALPGEN. The measured 
cross section exceeds the NLO prediction by about three 
o~ standard deviation. 



12. Conclusions 

Many measurements of vector boson plus light and 
heavy flavour jet production have been presented. Per- 
turbative QCD predictions by means of MCFM [JJ are 
in good agreement with data in all inclusive vector bo- 
son plus jet production measurements. The predictions 



of ALPGEN and MCFM are about three cr below the 
exclusive W boson plus b jet cross section measurement 
of CDF, where a reliable secondary vertex tagger has 
been used. The jet over Z boson pj balance measure- 
ment of CDF shows up limitations for the Jet Energy 
Scale precision at LHC experiments of the order of 3%. 
The pj imbalance can be explained by large angle fi- 
nal state radiation in form of collinear sub-leading jets. 
There is no perfect Monte-Carlo event generator. This 
holds for PYTHIA and HERWIG+JIMMY as well as 
for SHERPA and ALPGEN which are superior to the 
former parton shower Monte-Carlo event generators. In 
general the data are corrected for detector effects, so that 
they can be compared to predictions at the hadronic fi- 
nal state. This means that the data can be re-used for the 
tuning of Monte-Carlo event generators at any time in 
the future. 
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